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ABSTRACT 


A summary of modifications to the JANNAF version of Aerotherm's Boundary Layer Inte- 
gral Matrix Procedure (BLIMP-J) is presented. The BLIMP-J program is the standard prediction 
method for boundary layer effects in liquid rocket engine thrust chambers. One of these modi- 
fications, the addition of the turbulent models of Cebeci and Smith and of Bushnell and Beck- 
with, is discussed in detail. Comparisons of the predictions of these models and the model 
of Kendall between the models and with data are also presented. The Cebeci-Smith and Bushnell 
Beckwith models are in close agreement for all cases studied. For the case of large values 
of T g /T w (=5) the Kendall model does not agree with the other two. It is recommended that 
the Kendall model remain as the turbulent model for use in the standard prediction procedure 
until sufficient comparisons with appropriate high temperature ratio data are made. Recom- 
mendations for future work are alos presented. 
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SECTION 1 


INTRODUCTION 

Accurate predictions of the thrust loss due to boundary layer effects and of the wall 
heat flux are very important to the design and performance evaluation of rocket nozzles. The 
JANNAF version of Aerotherm's Boundary Layer Integral Matrix Procedure (BLIMP-J) is the reference 
method in the recommended rigorous analytical procedure for rocket engine performance predic- 
tion and evaluation (Reference 1). The JANNAF version of the BLIMP program has been specially 
modified for easy interface with other programs of the JANNAF recommended procedure. The 
BLIMP-J program is a fast and accurate procedure for solving the set of boundary layer equa- 
tion (momentum, energy, and species) for laminar or turbulent, chemically reacting flows with 
a wide variety of boundary conditions. This computer program has been fully documented and 
its use described in Reference 2. 

In this report the most recent improvements to the BLIMP-J program are described and 
the results of comparison of the turbulent models of Kendall (References 3 and 4), Cebeci and 
Smith (References 5 and 6) and Bushnell and Beckwith (References 7 and 8) are presented. 

Section 2 contains a summary of the recent modifications and additions to the program. The 
turbulent models and the results of the comparisons are given in Section 3. Conclusions and 
recommendations for future work are presented in Sections 4 and 5. 




SECTION 2 


CODE MODIFICATIONS 

This section contains a brief discussion of the additions and modifications made to 
the BLIMP-J program since the conclusion of a previous updating (Reference 9). Details of 
these modifications and the entire BLIMP-J program can be found in Reference 2. 

2.1 HOMOGENEOUS GAS OPTION 

It is frequently desired to make fast preliminary boundary layer calculations or to 
make calculations for air flow models where no chemical reactions take place. To facilitate 
these calculations a homogeneous, or nonreacting, gas option has been added to the program. 

In this option the gas species composition is specified by the user and is not allowed to 
change. Additionally, the gas viscosity and Prandtl number as functions of temperature are 
specified by the user. In this way the species equations are eliminated from the set of boun- 
dary layer equations and the complex equilibrium chemistry and transport property routines are 
replaced by very simple calculational procedures. This can greatly reduce the execution time 
and, in the case of low temperature air, provide greater accuracy through improved transport 
property values (p, Pr) at low temperatures. 

2.2 BINARY DIFFUSION OPTION 

Another simplification useful for reducing the execution time for problems with more 
than two elements in the chemical system is the binary diffusion option (see Section 2.4 and 
6.4 of Reference 2). The use of this option reduces the number of species equations to one 
from a maximum possible of seven. This can result in a 40 to 75 percent savings in total 
execution time, depending on the number of elements in the gas. (A typical problem might 
require 20 to 30 minutes (Univac 1108} without the binary diffusion option.) 

The effect of this option on the solution depends on the degree to which unequal mo- 
lecular diffusion is important to the accuracy of the solution. For turbulent flow problems 
with no blowing the effects may be, and usually are, small. 
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2.3 LOW TEMPERATURE RANGE 


In certain cases where very large expansion ratios are encountered, it is sometimes 
the case that the gas temperature may vary over a wide range. A third range for thermody- 
namic curve fit constants has been added to allow for more accurate thermodynamic properties 
determination. This range is only used when specifically requested by the user. 

2.4 EXPANDED NAMELIST INPUT 

The input for the BLIMP program has been modified to accept namelist input for most of 
the input. This option facilitates minor changes to data decks as an aid to the user. Many 
input values have been preprogramed for standard default values when used with the namelist 
options. The complete formatted input has also been retained. 

2.5 MISCELLANEOUS MODIFICATIONS 

2.5.1 Plot Output 

Optional output to a plot file of many key parameters has been added. The boundary 
layer edge values, integral parameters, and boundary layer profiles can be stored (either 
on the file or on tape) for use with a plot routine. No plot routine is supplied with the 
program. The exact format of the plot file is described in Section 7.5 of Reference 2. 

2.5.2 Boundary Layer Flow Rate 

The mass flux in the boundary layer has been added as output. This was done so that 
comparison could be made with the mass flux in the wall zone of the TDK (Reference 10) solu- 
tion. It is desirable to have the edge of the boundary layer contained within the wall zone 
of the TDK solution. 

2.5.3 Edge Conditions 

The edge conditions for the boundary layer solution were previously calculated from a 
given stagnation condition and edge pressure distribution for an isentropic expansion. This 
has now been expanded to allow for a nonisentropic expansion whereby the edge thermodynamic 
state can be calculated from either pressure and velocity or pressure and entropy. In addi- 
tion the edge pressure gradient and edge velocity gradient can be input directly if desired. 

2.5.4 Convergence Criteria 

The convergence criteria for the boundary layer equations was reformulated to be self- 
adjusting for boundary layers with large streamwise pressure gradients or rapidly changing 


4 



boundary layer shape and thickness. These criteria compare the errors in the momentum, 
energy, and species equations (see Section 3 of Reference 2} to a relatively small number 
which is adjusted according to certain parameters in the equations. The present methods have 
proved to be very reliable for all types of problems considered. 
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SECTION 3 


' TURBULENT MODELS STUDY 

The prediction of wall heat flux, momentum thickness, and displacement thickness for 
turbulent boundary layers is very dependent upon the turbulent model employed in the predic- 
tion tool. The original turbulent model in the BLIMP code was developed by R. M. Kendall, 
has been described in previous BLIMP documentation (References 4 and 11), and has extensive 
use for both external and internal flow predictions. There are many other models in use for 
turbulent predictions. Two of these, by Cebeci and Smith (References 5 and 6) and Bushnell 
and Beckwith (References 7 and 8), were selected to be added to the BLIMP code. This would 
make possible a comparison of three of the leading turbulent models of the mixing length type 
in the same computational routine for general equilibrium chemistry and a broad range of pos- 
sible boundary conditions. 

In Section 3.1, the basic formulation of each of the three turbulent models is given. 
In Section 3.2, a comparison of the three models is presented. The execution time for the 
three models is compared in Section 3.3. 

3.1 THE TURBULENT MODELS 

The general features of the turbulent models are described in Section 3.1.1 and each 
model is described in detail in the subsequent sections. 

3.1.1 General Features 

Bossinesq's eddy viscosity concept is adapted to write the Reynolds stresses as 


- (pv)’u' 


pe 


jju 

m 9y 


and a similar relation is used to define eddy conductivity, e h 


( 1 ) 


- (pv)'T' = 


pe. 


ST 
h 3y 


( 2 ) 


PRECEDING PAGE 


BLANK NOT FILMED 


7 



All the three models under present discussion employ the Prandtl mixing length hypo- 


thesis in which it is assumed that 


= £v , 


(3) 


where l is the mixing length and v^ is the turbulent velocity. The differences between the 
three models comes about through the formulation of l and v^. Kendall and Cebeci treat the 
boundary layer as a composite layer consisting of inner and outer regions. In the inner re- 
gion the turbulent velocity is written as 


v 


t 



(4) 


and the mixing length is assumed to be proportional to the distance from the wall. In the 
outer or wake region, the boundary layer is assumed to behave similarly to free turbulent 
shear flow with v^. = u g , the free stream velocity, and i l = c6* where c is a constant and 6* 
is a boundary layer characteristic thickness taken as the velocity defect thickness. Thus, 


where 


t) ay 


(6) 


Bushnell and Beckwith, however, treat the boundary layer as a single layer and apply Equation 
(4) throughout by introducing the intermittency concept in the definition of H. The most fun- 
damental differences in the models arise from the manner is which the mixing length expression 
is obtained. The Cebeci and Bushnell expressions originate from Prandtl 's proposal that in 
the region of the development of turbulence. 


dS, _ 
dy 


k 


(7) 


which has solution 


l = ky 


( 8 ) 
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From this solution many significant modifications have followed to account for the effect of 
variable properties, pressure gradient, Reynolds number, etc. It is significant that these 
modifications were made to the solution and not to the basic proposition as expressed by Equa- 
tion (7). The Kendall model, on the other hand, originates from modifications to the basic 
proposition to account for the effects of variable properties (see Equation (10)). Thus, dif- 
erences in the models might become more pronounced as the degree of property variation in- 
creased. This is exactly what is reflected in the comparisons presented in Section 3.2. 

The turbulent transport of scalar quantities is treated the same way as momentum by 
introducing the concepts of eddy conductivity, e^, and eddy diffusivity, e^. Turbulent Prandtl 
and Schmidt numbers are defined as Pr^ 5 e m /e^ anci Sc t 5 e n/ e D' Ce b ec ^ Proposes an expression 
for Pr^ as a function of the distance from the wall but in the Kendall and Bushnell -Beckwith 
models Pr^ is assumed to be a constant. The turbulent Schmidt number is also taken to be con- 
stant in all the models. 

3.1.2 Kendall Model 

This model employs the two-layer concept of the turbulent boundary layer. The wall law 
is based on the following three concepts: 

• lim l = 0 
y -*■ o 

• lim dJl/dy = 0 
y + o 

• Rate of increase of the mixing length with y is proportional to the difference be- 
tween the value postulated by Prandtl (ky) and its actual value 

The proportionality factor in this relation is assumed to be dependent on the local shear 
stress and local kinematic viscosity 

^=(ky-*)^R (9) 

y V 

where y + is a constant. The values of the constants k and y + recommended in this model are 
0.44 and 11.823, respectively. These constants have been obtained by matching the predictions 
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with experimental data in incompressible turbulent boundary layers with and without blowing (Re- 
ference 11). (Physically k can be considered as a measure of the rate of growth of the mixing 
length with respect to distance from the wall and y* is a measure of the thickness of the lami- 

a 

nar sublayer.) The validity of the model for flows with wall blowing and streamwise pressure 
gradient is argued on the basis of using the local flow properties, such as local shear, in 
the model. 

For compressible flow, the wall law is modified as follows: 


dpt _ 

dy 


[k 



pdy - pi.] 


/E7p 


( 10 ) 


where, instead of describing the length scale of a turbulent eddy, the mass of the eddy, pt, 
is related to the mass available between the wall and the point of interest. The constants 
k and y + , however, are left at their incompressible values. The above integral -differential 
equation is solved numerically to obtain the local value of the mixing length £. 

In the wake region, it is assumed that the eddy viscosity is a constant and is given 

by Clauser's expression (Equation (5)) where c = 0.018. The wall and the wake regions are 

matched by the following procedure: the expression for the wall region is used until it 

exceeds the wake value at which point the wake value of e is used for the remainder of the 

m 

boundary layer thickness. This value is linearly damped in the outer-portions of the boundary 
layer so that a value of zero occurs at the boundary layer edge. 

3.1.3 Cebeci -Smith Model 

As was mentioned above, a two-layer model is also used by Cebeci and Smith. In the 
inner (wall region) the Van Driest (Reference 12) form of mixing length is now used: 

£ = k^ [1 - exp(-y + /A + )] (11) 


where 

y + = y A w /p /v 

Van Driest suggested constant values of 0.4 and 26 for the k m and A + , respectively. (These 
have essentially the same meaning as k and y* and were arrived at by comparison with incom- 

a 

pressible air flow data.) In the Cebeci model, however, these constants are replaced by 
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functions accounting for pressure gradient, and blowing. Compressibility effects are also 
accounted for by using local values for p and p. 

For flows with pressure gradient and mass transfer, Cebeci replaced the wall shear in 
the damping parameter by x s which he obtained from the simplified form of the momentum equa- 
tion in the sublayer (Reference 5): 


= df> 

dy v w s dx 


The solution of this equation at y =11.8 results in 


A + = A ^ [exp(ll .8 V*) - 1] + exp(ll .8 V+l 
^w 


(12) 


(13) 


where 


+ 

P = 


dx 


P u 3 

P W T 



(14) 


and A = 26. 

Following Van Driest 's approach to arrive at the mixing length formulation with a 
damping factor in the inner region, Cebeci derived the following expression for eddy conduc- 
tivity: 

e h = k m k h y 2 [1 " exp(-y + /A + )][l - exp(-y + fPr/B + )] 

where 

B + = tsj^r [exp(i 1 .8 V*) - 1] + exp(11.8 V*)j“ V2 (16) 

V w 


Ju 

3y 


(15) 


and 


k h = 0.44 , B = 34 

Cebeci (Reference 6) further agreed that the above values of k^, k^. A, and B are only 

satisfactory for large Reynolds number (Re„ > 6000), and he proposed functions of Re. for k , 

u urn 

k , A, and B. There is some question as to the validity of the Re. dependence, particularly 

h . ■ “ 

*If V + = 0, A + = A {-11.8p + + 1}' 1/2 
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for compressible flows in nozzles. Furthermore, the model is completely adequate without this 
dependence;* for these reasons the constant values of k m , k^, A, and B are used.** 

The turbulent Prandtl number (Pr t = e m /e h ) is obtained from Equations (3), (4), (11), 
and (15); 

k [1 - exp(-y + /A + )] 

p r = 0! C171 

1 k h [1 -exp(-y + /PF/B + )] ( ) 

Although Equations (11) and (15) are valid only in the boundary layer inner region, Cebeci 
shows that Equation (17) agrees satisfactorily with experimental data of References (13) 
through (15) throughout the boundary layer and, hence, it is so used. 

In the wake region, Cebeci (Reference 6) uses the Clauser expression for eddy viscos- 
ity, Equation (5) with c = 0.0168. This expression is damped in the same way as in the Ken- 
dall model. 

3.1.4 Busline! 1 -Beckwith Model 

The Bushnell-Beckwith model is a single layer model which reduces to the Van Driest 
form of mixing length near the wall and is modified in the outer region by an intermittency 
factor y (Reference 16). The mixing length expression is written as: 

| = K [1 - exp(-y + /A + )] f(y/6b Vi! (18) 


where 


and 


- 1 - erf [5(y/<S - 0,78)] 
2 


(19) 


f(y/6) = tanh(jl^) (20) 

and the constants are: k^ = 0.4, K = 0.08, A + = 26 (from incompressible flow comparisons). 

The boundary layer thickness 6 appearing in Equations (18) through (20) is defined as the 
distance normal to the wall where the velocity ratio (u/u g ) = 0.995. 

"k , . 

Personal communication: Tuncer Cebeci, McDonnell -Douglas, Long Beach, California. 

**In a later article (AIAA, Vol. 11, No. 1) Cebeci argues that k m and A are constant. 
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The present model has been tested against experimental data by Bushnell and Beckwith 
(Reference 17). In their work, however, they use a different function than the one given by 
Equation (20). They assume that f = y/6 in the inner wall region, y/6 <l 0.1, and it is a func- 
tion of the incompressible shape factor (H = 6*/e) in the far wall region, y/6 > 0.3. The 
values of f in the far wall region are obtained from a curve fit to experimental data of Sl/6 
versus H. In the interval 0.1 < y/6 < 0.3 a straight line is used to join the inner and far 
wall regions. Based on this model, Bushnell and Beckwith compared their predictions of flows 
with blowing and pressure gradient with experimental data of References 18 and 19. They re- 
port (Reference 17) that in the application of the model to flows with wall blowing, the ef- 
fect of blowing could be accounted for only when the wall damping factor of Van Driest, A + , 

was made an experimentally based function of the blowing rate. The present functional form 
of f, Equation (20), is based on the recommendation of Harris.* 

As it is noted by Harris (Reference 8), based on the available data, there exists a 
lack of conclusiveness as to how the turbulent Prandtl number should be formulated in terms 
of local boundary layer parameters under different flow conditions. Therefore, a constant 
value of 0.9 is used for Pr t in this model. 

3.2 COMPARISONS OF THE TURBULENT MODELS 

The BLIMP program has been previously verified for several sets of turbulent flow data 

(Reference 9). In these cases the turbulent model of Kendall was used. The predictions were 

in good agreement with the data for the air flow cases. Two of the cases reported in Refer- 
ence 9 are used here for comparison of the turbulent models. In addition, several hypotheti- 
cal cases are evaluated to examine potential differences in the models and one set of data for 

flow over a hemisphere cylinder with large values of T /T (»5) is discussed 

w w 

3.2.1 Back and Cuffel — Supersonic Nozzle with Heat Transfer 

These experiments were carried out in a cooled, conical nozzle with a convergent and 
a divergent half angle of 10° at the Jet Propulsion Laboratory, California Institute of Tech- 
nology (Reference 20). The air flow was tripped well upstream of the converging section so 
that the flow was fully turbulent throughout the nozzle. These flow conditions are relevant 
to conditions in rocket engines operating at thrust levels for which laminarization does not 
occur. Values of T g /T w vary between 0.7 and 2; therefore, density variations through the 
boundary layer should not cause significant differences between the models. Values of Re. 
were between 6000 and 15,000 and the Mach number at the exit was 3.55. 

* 

Personal communication. 
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Static pressures, wall heat fluxes, and coolant-side wall temperatures were measured along 
the nozzle wall. Heat flux was determined by calorimetric measurements in circumferential cool- 
ant passages. Boundary layer surveys were made with a flattened pitot probe 0.000127 m (0.005 
inch) in height and with thermocouple probes at five locations, two of which are shown in Fig- 
ure 1. 

Comparison between the mode's and, where possible, the data for the two stations shown 
in Figure 1 and for the wall heat flux throughout the nozzle are presented in the following 
paragraphs. All of the model predictions were done with identical conditions using the re- 
commended values for the various input parameters. In addition, the Cebeci model was also 
used with variable turbulent Prandtl number. 

Figure 2 shows the overall agreement of the models and the data as reflected by the 
wall heat flux. The only significant difference is for the variable turbulent Prandtl number 
case in the converging section of the nozzle. In Figure 3, the momentum thickness is shown 
for the various models. Again, the overall agreement is very good. Tables 1 and 2 and 
Figures 4 and 5 show the comparisons near the entrance to the converging section and near the 
exit to the diverging section. The agreement among the models and the data is again very good. 

In almost all instances the models are within 10 percent of each other, with some exceptions 
for the Cebeci model with variable turbulent Prandtl number. The heat flux to the wall is an 
indication of the near wall behavior and the momentum thickness reflects the general profile 
shape. Thus, the models appear to be in agreement in both respects. This is further evi- 
denced by comparison of the eddy viscosity profile shown in Figure 6. Here it can be seen 
that throughout the portion of the boundary layer dominated by turbulence (c > v) the models 
are in very close agreement. The differences in the laminar sublayer (v > e) have little ef- 
fect on the solution since the molecular viscosity dominates. The slight differences in the 
outer wake region also have little effect since the shear terms are small as a result of 
du/dy being very small. 

3.2.2 Rocketdyne: Og/Hg Two-Dimensional Nozzle 

This case* is representative of the type of data that can be expected from hot fired 
nozzles using Og/H^ fuel systems. The gas side wall temperature and the wall heat flux dis- 
tributions were measured; however, no boundary layer measurements were made. Since heat flux 
is an important quantity and this case was for a representative liquid rocket fuel, it was 

* 

The data for this case were supplied by Mr. George Osugi of Rocketdyne. 
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Figure 3. Comparison of momentum thickness for Back and Cuffel data 
R T = 0.0202 m (0.06625 ft) 



TABLE 1. BACK AND CUFFEL CASE, x/R 
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Figure 5. Comparison of velocity profiles for Back and Cuffel data at x/Rj = 10.797 

u e = 1100.9 m/s (3612 ft/s) 

f) = 9.652 x 10' 3 m (3.167 x 10~ 3 ft) 









felt that comparative predictions would be of interest. Also, density variations (caused by 
temperature and species composition changes) of about a factor of 4 occur across the boundary 
layer and this was expected to be a possible cause for difference between the models. Values 
of Re. were between 1300 and 3600, much lower than the 4000 - 20,000 expected in SSME. It 
should be noted that this data is subject to uncertainty in the mixture ratio which make ac- 
curate comparison with predictions very questionable. 

The nozzle geometry, fuel mixture ratio, stagnation conditions, calculated axial pres- 
sure variation, and wall temperature variation were provided. (The nozzle is 0.0762 m (0.25 
ft) wide.) The stagnation conditions are given as: 

P Q = 4.6182 x 10 6 N/m 2 (45.57 atm) 

T q = 3570°K (6430°R) 

MR = 6.15 

The injector plane was at the entrance to the nozzle and the initially low heat trans- 
fer shown by the data (Figure 7) was assumed to result from the presence of a liquid layer 
near the injector. Accordingly, the prediction was started downstream of the injector and 
assumed to have an established boundary layer at the starting position. The mixture ratio 
in the boundary layer was assumed to be 6.15. The heat flux data and the results of the 
BLIMP predictions are shown in Figure 7. The Kendall model is about 15 percent higher than 
the Bushnell model and the Cebeci model is about 6 percent lower. At the exit Kendall is 
about 38 percent higher and Cebeci is about 4 percent lower than the Bushnell. The predic- 
tions are for a fully turbulent boundary layer with a starting length of 0.0561 m (0.18406 ft). 
It should be noted that the flow was suspected of being of the laminarizing type. Also, the 
predicted values are sensitive to the mixture ratio. If the wall zone is fuel rich (lower 
mixture ratio) the predictions are uniformly lowered. The starting length is also an influ- 
ential parameter. For these reasons no attempt was made to achieve better agreement with the 
data. Instead, the case was used to compare the three turbulent models for this type of flow, 
e.g., large edge temperature to wall temperature ratio (5 to 6) and chemical reactions. It 
is quite clear that the three models do not agree for this case. The velocity profiles are 
shown in Figure 8. The Kendall model is the most different of the three. This is attributed 
to the fact that the wall laws for the Bushnell and Cebeci models are very similar (of the 
Van Driest type) whereas the Kendall model is basically different, consisting of an integral- 
differential formulation. As discussed previously, Section 3.1.1, it was expected that these 
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Bushnel 1 



1300 - 3200 





differences might become significant for large property variations through the boundary layer. 
The differences in the models are also apparent in Figure 9, which shows the eddy viscosity 
profiles. (The degree of property variation across the boundary layer is also indicated by 
the edge and wall values of kinematic viscosity.) 

3.2.3 Hypothetical Cases 

As a consequence of the differences observed in Section 3.2.2, the models were compared 
for several cases which would help to define the differences in the models. These cases and 
the results will be briefly summarized. 

The comparison in Section 3.2.2 was for relatively low values or Re 0 and a severe pres- 
sure gradient. Therefore, the SSME flow conditions were used to compare the models at high 
Reynolds number (2000 < Re Q < 20,000) and less severe pressure gradient. Large values of 
T fi /T w were used. The results were in agreement with those of Section 3.2.2. The Kendall 
model was between 25 and 40 percent higher in wall heat flux than the other two models. 

The effects of property variation only were examined by using a zero pressure gradient, 

Re„ = 6000 air flow case for T /T =5 and T /T =0.2. For the large values of T /T the re- 
0 ewew ew 

suits were again similar to the Rocketdyne and SSME cases. However, for T./T “ 0.2 the models 
were in very close agreement. Examination of the profiles (for T /T ~ 0.2) revealed that most 

6 W 

of the property variation occurred in the wake portion of the boundary layer. (Viscous heating 
maintained the temperature near the wall value in the inner region of the boundary layer.) 

Since the Kendall and Cebeci models are identical in the wake no differences should occur if 
the wall region is approximately a constant property region. This further confirmed that it 
is for large property variations in the law of the wall region that the models disagree. 

As a further test for large property variations with large Re 0 , the Back and Cuffel 
case was rerun with the wall temperature reduced so that T g /T w = 5. The results were the 
same as in all the previous cases of large T g /T w . 

3.2.4 Flow Over a Hemisphere-Cylinder 

At the time the turbulent models were being examined some experimental data from AEDC* 
were made available. This data was for M = 12 air flow over a hemisphere (0.0635 m radius) 
with a cylindrical afterbody. Over the region of interest Re. was less than 300. The body 

Unpublished data, provided courtesy Mr. Emmitt Edenfield, Arnold Engineering and Development 
Center. 
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shape and pressure distribution are shown in Figure 10. The interesting feature of this 

data was that the ratio T /T varied from 5.2 to 2.0. BLIMP predictions for this data using 

e w 

an external flow option in another version of BLIMP were made. The flow over the hemisphere 
was artificially tripped at S/R Q * 0.2 on one side of the hemisphere. (The BLIMP predictions 
for heat flux for the laminar flow along the untripped side were in excellent agreement with 
the data.) The heat flux predictions for the turbulent flow for the Kendall and Bushnell 
models are shown in Figure 11. The Cebeci model is approximately the same as the Bushnell 
model. In the region of large T g /T w the Kendall model is as much as 40 percent larger than 
the Bushnell model. 

This single comparison is not considered as conclusive because the data is not of the 
type desired for verification of a turbulent model. First, the models are not considered 
to be well tuned for low Reynolds number flows and second, there is some evidence that in the 
region immediately downstream of transition there is an increase in turbulent mixing that has 
not been taken into account in the models. Although part of the differences in the predictions 
of the models can be attributed to the severe property variations, the agreement with the data 
by the Kendall model may be fortuitous. 

3.3 COMPARISON OF COMPUTATION TIME 

There are many factors that effect the execution time of each model - the starting 
profile, the degree of nonsimilarity, etc. Each model performs differently in any given 
situation. For the several problems considered there was no clear indication of which model 
is fastest. The time per iteration is usually slightly more for the Bushnell model (~10 per- 
cent) than for the Cebeci or Kendall models, which are about the same. However, the number 
of iterations does not consistently favor any of the models. For example, in the three cases 
initially considered (JPL data, Rocketdyne Og/H., case, and SSME model air flow) a different 
model had the shortest overall execution time for each case. Thus, there does not appear to 
be a strong motive to select any model based on computation time. 


27 



Figure 10. Body shape and pressure distribution for AEDC hemisphere 
cylinder. 

M = 12 

R 0 = 0.0635 m (2.5 in.) 

P 0 = 1.94 x 10 5 N/M 2 (1.912 atm) 







lux for AEDC hemisphere cylinder. 

I 

.0635 m (2.5 in.) 

= 8.74 x 10 5 watts/m z (77 Btu/ft 2 s) 
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SECTION 4 


CONCLUSIONS 

The major conclusions of this work pertain to the turbulent model study and are given 

below. 

• For nonreacting flows with moderate pressure gradients (i.e., not severe enough 
to cause laminarization) and values of T g /T w < 3 the three models are equivalent 
and seem to agree with the data. 

• For nonreacting flows with large values of T g /T w (“5), the Kendall model does not 
agree with the Cebeci or Bushnell models. The Kendall model predicts significantly 
higher values of heat flux. 

• The Kendall model should remain as the standard model for prediction of boundary 
layer effects in liquid rocket engine thrust chambers since it is the model that was 
initially incorporated in the BLIMP program. It is imperative that data be obtained 
to resolve the question of large property variations in the "law of the wall" region. 
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SECTION 5 


RECOMMENDATIONS 

1. The importance of gas phase kinetics on the prediction of wall heat flux thrust loss, and 
body displacement should be assessed. If kinetics prove to have significant effects, a 
means to include kinetics in BLIMP-J should be investigated. 

2. The capability to freeze the chemical composition at any point in the expansion should be 
added to BLIMP-J. Currently BLIMP-J can be run either with a completely equilibrium 
chemistry or with frozen chemistry. This option would allow the user to mix the two op- 
tions. For example, equilibrium chemistry might be used upstream of the throat, and the 
composition frozen at the throat values for the expansion. 

3. A parametric study to determine the sensitivity of the predictions to certain input values 
should be made. Two of the most obvious parameters are the mixture ratio and the initial 
starting length for the boundary layer. 

4. The importance of laminarization for the rocket nozzles of interest should be assessed. 

The current turbulent models do not appear to adequately model this effect; although a 
rigorous comparison with data has not been made. If laminarization is important, adequate 
means to model it should be investigated. This could include modification to the current 
models in BLIMP or the addition of a suitable turbulent kinetic energy model. 

5. In view of the differences in predictions of the three turbulent models for the case of 
large T e /T w> suitable data should be identified and comparisons made. If no such data are 
available then an experimental program should be designed and conducted. 
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